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We present a scenario for the peculiar coexistence of charge fluctuations observed in quasi- 2D 
1/4-filled organic conductors #-(BEDT-TTF)2X in the quantum critical regime where the charge 
ordering is suppressed down to zero temperature. The scenario is explored in the extended Hubbard 
model including electron-phonon couplings on an anisotropic triangular lattice. We find that the 
coexisting fluctuations emerge from two different instabilities, the "Wigner crystallization on lattice" 
driven by the off-site Coulomb repulsion and the charge- density- wave formation due to the nesting 
of the Fermi surface, not from phase competition or real-space inhomogeneity. This mechanism 
explains the contrastive temperature dependence of two fluctuations in experiments. 
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Charge ordering (CO) is a periodic arrangement of 
electrons which drives a metal-insulator transition. It 
is often found in organic conductors as well as transition 
metal compounds, and it plays a key role in electronic 
properties such as transport phenomena. Careful control 
of chemical substitutions and/or the external pressure 
has revealed a systematic change of the electronic states, 
and has opened up the possibility of a comprehensive 
understanding of the correlated electron systems p], Q • 

Among them, quasi-2D organic conductors 6-(ET)2X 
are intriguing systems, where ET is an abbreviation for 
BEDT-TTF and X represents monovalent closed-shell 
unit. The crystal structure consists of an alternating 
stack of ET and X layers; ET molecules constitute con- 
ducting layers of an anisotropic triangular structure as 
schematically shown in Fig.[TJ and anions X~ form insu- 
lating layers. The charge transfer between ET and X~ 
layers leads to hole doping at 1 /4 filling in the ET layers 
(one hole per two ET molecules on average). By chang- 
ing X units and/or by applying pressure, the compounds 
exhibit a rich variety of electronic phases at low temper- 
atures (T), such as CO insulator, spin-Peierls state, and 
superconductivity. It was proposed that the variety is 
summarized in a systematic way by using the angle be- 
tween neighboring ET molecules, i.e., the dihedral angle 
(j) (Fig. HJ) |3, |j|. It is considered that <p modifies the 
transfer integrals between ET molecules, t p and t c , in 
particular, the ratio between them, t c /t p . In fact, in the 
series of X = MM / (SCN) 4 with M' = Co or Zn, a stripe- 
type CO takes place in compounds with M = Tl and Rb 
H, @, 0, @1 , and the transition temperature is decreased 
from T co = 250 K for Tl to 190 K for Rb, as <j> de- 
creases from (j) = 116° to 111° and correspondingly t c /t p 
decreases from 0.5 to 0.4. Finally, Tco is suppressed 
down to the lowest T for M = Cs with <\> = 104° and 
t c /t p ~ 0.1. 

The marginal compound with M = Cs, close to the 
quantum critical point where Tco — ► 0, has attracted 
special attention because of its peculiar low-T proper- 
ties. A striking feature was found in X-ray scattering re- 
sults: two diffuse peaks were observed at different wave 



numbers 0, [TO]. One appears at the same wave num- 
ber as that of the stripe- type CO in Tl and Rb com- 
pounds, which has a twofold period along the c axis. The 
other peak corresponds to a threefold period both in the 
a and c axes. It is remarkable that charge fluctuations 
with different wave numbers coexist. Moreover, these dif- 
fuse peaks exhibit contrastive T dependence: The peak 
for the twofold period grows rapidly below T ~ 90 K, 
while the threefold one shows little T dependence [0, E3] • 
Recently, anomalous transport properties were also ob- 
served [ll|, ll2|, ll3[ , and their relation to the characteristic 
charge fluctuations has been issued. In fact, it was found 
that the peak for the twofold period is suppressed under 
the electric current, while the threefold one is not [l2[. 

A scenario to understand the peculiar behaviors is 
phase competition among different CO states. Several 
different CO states are stabilized due to strong Coulomb 
repulsion [14], ll5|, [lg, [l7| , which are sometimes considered 
as generalization of "Wigner crystals" [18]. It was shown 
that, in a certain parameter regime, some of the CO 
states are energetically degenerate, and a keen competi- 
tion among them may su ppr ess long-range CO and even 
lead to metallic behavior jl9[ I20ll21, 22, 23]. It was spec- 
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FIG. 1: Schematic picture of BEDT-TTF layer. Ellipses rep- 
resent BEDT-TTF molecules, which form an anisotropic tri- 
angular lattice. (j> is the dihedral angle, and t p and t c are the 
transfer integrals for the nearest-neighbor molecules used in 
the tight-binding Hamiltonian. See the text for details. 



2 



ulated that the phase competition leads to an inhomoge- 
neous state where the competing phases are mixed in the 
real space by forming domainlike structure jll|, [l2[ . How- 
ever, it is not obvious how such inhomogeneity emerges in 
the clean system almost free from impurities and why dif- 
ferent domains grow with the contrastive T dependence. 
Furthermore, there is no experimental evidence for the 
domainlike state so far. 

In this letter, we present an alternative scenario for 
the puzzling coexistence of charge fluctuations as well 
as the systematic change of the CO phase diagram in 
#-(ET)2X by focusing on different roles of the off-site 
Coulomb repulsion and the kinetic energy of electrons. 

Our Hamiltonian consists of two contributions as H = 
Hehm +7^ei-ph- The former term represents the extended 
Hubbard model in the form 

hehm = t p ( c L c 3* + H - c -) + ^ ( c t c 3* + H - c ) 

(i,j)p,<r (hj)o<r 

+ U ^^n^riii + V ^2 n % n r (1) 

Here, the summations with (i,j) p ( c ) are taken over 
the nearest-neighbor sites along the x, y axes (c axis) 
(Fig. [1]). The off-site Coulomb interaction V is taken 
to be isotropic because of the almost equal distances be- 
tween neighboring ET molecules. We consider the case 
of t p , t c > and the average electron density at 1 /4 fill- 
ing, which corresponds to 6-(ET)2X by the electron- hole 
transformation. We take t p = 1 as an energy unit, and 
set the Boltzmann constant /cb = 1. 

The latter term 7Y e i- p h describes the electron-phonon 
couplings which are necessary because CO appears with 
a change of the lattice structure 0, 0, [Hj]. Although 
realistic phonon modes are complicated, here we sim- 
ply incorporate the Holstein type and the Su-Shrieffer- 
Heeger (SSH) type phonons in the model as H e \- P h = 

^Holstein + ^SSH with 

^Holstein = 9^ U ^ + E { ^jT U i + } ' ^ 

i i 

HsSH = 9p E ^ Ua ^ ~ U a,i)( c L C j° + Hx 
+ 9c ^ ~ U ^i)( C \a C jcr + H.C.) 

+ £ {^«« + <*) + + }• (3) 

i 

In 7^Hoistein, 9 is the coupling constant; U{ and de- 
note the lattice distortion and its conjugate momentum 
at site z, respectively. In 7Yssh, we describe the displace- 
ments of zth site in the a and c directions by u a ,i and 
u Ci i, respectively, and we assume that they modulate the 
transfer integrals t p and t c independently with the cou- 
pling constants g p and g c \ the summations are taken in a 
manner that an elongation (shrinkage) of the bond length 



decreases (increases) the transfer integrals. U a i and H c i 
are the conjugate momenta of u a ^ and u c ^ respectively. 
M, M' and c^o, u; f are mass and eigen- frequency of the 
corresponding phonons, respectively. 

To study CO phenomena in our model, we mainly an- 
alyze the charge susceptibility Xc(q, iwq) by RPA. For 
instance, in the case of the model with the Holstein 
phonons (TYssh = 0), Xc is given in the form 

Xc (q, iwq) = : , ur , — : — \ — Tqw : — r , (4) 

where 

W(q, ioj q ) = U + 2V(q) + 2j 2 D^ (q, iw q ), (5) 
and is the bare susceptibility given by X^°^(q, ioo q ) = 

~k Ep{/(^P+q) - f(ep)}/{iVq ~ ( £ P+q ?p)h N is the 

number of sites, / is the Fermi distribution function, and 
e p — 2t p (cos p x + cos p y ) + 2t c cos(p x + p y ). Here, V(q) is 
the Fourier transform of the off-site Coulomb interaction, 
V(q) = 2V{cosq x + cos^ + cos(q x + q y )}; is the 

bare phonon Green's function given by (q, iu q ) = 
— 2uoq/(uo^ + cjq), and 7 = gj ' \[2MuJ§. The expression 
becomes more complicated when Hssu is included, but 
the main features of our results below are retained. 

Now we discuss RPA results, mainly for the model 
without the SSH-type coupling (7Yssh = 0) for simplicity. 
Figure[2](a) shows the phase diagram in the t c /t p -T plane 
for a typical set of parameters. Tqo is determined by 
the divergence of the static charge susceptibility, Xc(q, 0). 
The divergence appears at q = Q2 ~ (7r, 0) [and (0, tt)] 
although Q 2 weakly depends on t c /t p . Q2 corresponds 
to a stripe-type CO pattern with period two along the c 
axis as in the left inset of Fig. [U[b) (diagonal-type CO 
[15[). Tqo is lowered as decreasing t c /t p and goes to 
zero at t c /t p ~ 0.1, a parameter corresponding to the Cs 
compound in which Tqo — > is experimentally observed. 
Here, we note that the Holstein phonon stabilizes CO by 
effectively reducing the on-site repulsion U [251 ]. 

In the quantum critical regime at t c /t p ~ 0.1 where 
^co —> 0, Xc(q,0) shows two broad peaks at different 
wave numbers as shown in Fig. (2Kb): One is at Qi c± 
(27r/3, 27r/3) and the other is at Q2 ~ (tt, 0) [and (0,7r)]. 
Here the Qi peak corresponds to CO pattern with period 
three, as shown in the right inset of Fig. [51(b) [l7[. These 
two peaks exhibit contrastive T dependences as shown 
in Fig. (2](c): The Qi peak shows little T dependence 
in this T range, while the Q2 peak grows rapidly as T 
decreases. The behaviors are consistent with the X-ray 
scattering data in the Cs compound [9|, [l0[ . 

What is the origin of the coexistence of two charge fluc- 
tuations and their contrastive T dependence? The charge 
fluctuations at Qi and Q2 have different origins, the off- 
site Coulomb interaction V and the nesting property of 
the Fermi surface, respectively, as explained below. From 
eqs. and ([5]), it is clear that Xc(q 5 0) is enhanced at 
the wave numbers where the Fourier transform of the 
off-site Coulomb interaction V(q) is minimized or the 
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FIG. 2: (Color online), (a) RPA phase diagram for the model 
with the Holstein phonons at £7 = 3.2, V = 1.0, and 7 2 /ujq = 
1.5. (b) Charge susceptibility Xc(q 5 0) at t c /t p = 0.1 and 
T — 0.05. The insets show CO patterns corresponding to 
the peak wave numbers, (c) Temperature dependence of the 
charge susceptibility. The inset shows the peak heights, (d) 
Fourier transform of the off-site Coulomb interaction — V(q), 
and bare charge susceptibility x^(q, 0). The inset shows the 
parametrization of q = (q x ,q y ) in the plots of (c) and (d). 



bare susceptibility x^ (q? 0) is maximized. Figure [21(d) 
plots —V(q) and x^(q, 0). Obviously, the Qi peak cor- 
responds to the minimum of V(q), indicating that this 
peak originates from the off-site Coulomb repulsion. On 
the other hand, the Q2 peak corresponds to the max- 
imum of X^H^O), an d hence this originates from the 
noninteracting electronic state, i.e., the nesting property 
of the Fermi surface. These assignments explain the con- 



trastive T dependence of the two peaks: The Qi peak 
is of the Coulomb origin, and hence its height is weakly 
dependent on T in the present T range much lower than 
the energy scale of V. On the other hand, the Fermi de- 
generacy at low T <C t p enhances x^ (q? 0) at the nesting 
vector q = Q2, which leads to the rapid increase of the 
Q2 peak. In the inset of Fig.(2](c), we show T dependen- 
cies of the peak heights in a wider T range. The Q2 peak 
steeply grows at low T with a concave T dependence, 
while the Qi peak has a weak, convex T dependence and 
saturates at T ~ 0.1 [261 ]. 

Thus, our RPA results successfully reproduce the ex- 
perimental features in #-(ET)2X, i.e., the coexisting 
charge fluctuations and their contrastive T dependences 
as well as the systematic change of phase diagram in 
terms of t c /t p . The coexistence of charge fluctuations 
comes from the fact that the off-site Coulomb energy 
and the kinetic energy (the Fermi surface nesting) favor 
the different types of CO: the former induces an instabil- 
ity toward the "Wigner crystal" formation, whereas the 
latter leads to the charge-density-wave instability. The 
mechanism is intrinsic in the correlated electron systems, 
and it holds for the model including 7Yssh- 

There remain some issues to be clarified. One is the 
order of CO phase transition. The transition is always 
of second order in RPA, while it is of first order with 
a structural change in experiments 0, 0]. An interest- 
ing experimental observation is that a diffuse peak for 
the three-fold period appears above Xco [7, 8] with an 
unusual critical behavior in spite of the first-order transi- 
tion [5, 6]. A corresponding behavior is found in our RPA 
results; the threefold Qi peak develops above Xco since 
the peak is induced by V much larger energy scale than 
Xco- However, in order to describe the first-order nature 
of the transition, it is necessary to go beyond RPA. An- 
harmonic contributions from phonons may play a role as 
well. This is left for future study. 

Another issue is the wave numbers of CO and fluctu- 
ations. In our results, the coexisting charge fluctuations 
appear at Qi ~ (27r/3, 27r/3) and Q2 ~ (tt, 0). Whereas, 
in experiments, they are found at qi c± (7r, — 7r/3) and 
q 2 ~ (7r/2,7r/2) [SI, Glj. Accordingly, the stabilized 
CO pattern below Xco is also different; our result in 
Fig. [21(a) predicts the diagonal- type CO with Q2 ~ (7r, 0), 
while the experimental one is the horizontal type with 
q 2 ~ (7r/2,7r/2) 0, !, [3. We consider that these dis- 
crepancies are not fundamental and can be reconciled 
without modifying the essential part of our results in the 
following reasons. As to the Qi peak, the wave number 
is determined by the minimum of V(q), and it locates at 
q = Qi because we assume the simplest form of the off- 
site Coulomb interaction, i.e., the isotropic one for the 
nearest-neighbor sites only. If we include longer-range 
parts of and/or an anisotropy in the Coulomb interac- 
tion, V(q) takes its minimum at a different wave number, 
which may coincide with qi (2?J. On the other hand, the 
Q2 peak comes from the maximum of x^°\ an d hence, 
its wave number is dependent on the form of transfer 
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FIG. 3: (Color online). Phase diagram for the model with the 
SSH phonons at T = 0.01 and U = V = 7 = 0. CO patterns 
and displacements of ET molecules are also shown. 

integrals which modifies the Fermi surface nesting. In 
fact, we observe a shift of the peak when we include the 
SSH- type coupling which directly couples to the transfer 
integrals. 

To demonstrate how the SSH- type coupling affects the 
ordering vector Q2, we study the low-T phase diagram 
for the model including 7Yssh- We decouple Coulomb 
interactions and electron-phonon couplings by using the 
mean fields, (cl a c ia ), (c\ a Cj a ), (m) and (w a ( c ),;), and ob- 
tain a solution by minimizing the free energy among var- 
ious mean-field configurations within a unit cell of up to 
4x4 sites in the x and y coordinates. Figure [3] shows 
the result at T = 0.01 by changing the SSH couplings 
1p — dpi y/ZM'ujQ and j c = g c /^/2M f u f Q . To focus on 
the effect of TYssh, here we omit the Holstein coupling 



and take U = V = 0. A variety of phases emerges as 
shown in the figure; in particular, in a region where both 
7 P and j c are substantial, the horizontal-type CO is sta- 
bilized whose wave vector is consistent with q2 observed 
in Tl and Rb compounds [28]. We confirmed that the 
horizontal CO remains robust against finite U and V 
within the Hartree-Fock approximation. Thus, the SSH- 
type electron-phonon coupling plays an important role 
in determining an actual CO instability induced, which 
suggests possibility to shift Q2 to q2. For further com- 
parison, detailed information is indispensable on realistic 
electron-phonon couplings, which are possibly more com- 
plicated than our model. 

In summary, we have theoretically investigated the ori- 
gin of puzzling features in charge ordering and fluctua- 
tions in #-(BEDT-TTF)2X. The peculiar coexisting fluc- 
tuations are ascribed to two different instabilities, the 
"Wigner crystallization" driven by the Coulomb repul- 
sions and the charge- density- wave instability due to the 
Fermi surface nesting. Our scenario does not need to 
suppose either a competition among different "Wigner- 
crystal-type" CO in the strong-correlation regime nor 
real-space inhomogeneity. Our results provide a key to- 
ward understanding the low-T transport phenomena re- 
lated to the charge fluctuations. 
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